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*  Introduction  to  rockot  injectors 
»  Major  trends  since  the  end  of  the  cold  war 

deve,opments  in  rocket  injector  desii 

’  Case  stLJdy-'  gas/gas  injector  development 
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Introduction  to  rocket  injectors 


p  Major  kinds  of  rocket  injectors 
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Space  Shuttle  Main  Engine  Injector 
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From  Sutton,  “Rocket  Propulsion  Elements,”  e'"  ed.,  pg  277,  Wylie  1992 
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Platelette 


Face  plate  with  Fuel  distribution 

inclined  injection  holes  plate  with  etched 

(laser  drilled)  cooling  passages 


bonded 

injectors 


Cover  plate  with 
oxidizer  inlet 


From  Sutton,  “Rocket  Propulsion  Elements,"  6"’  ed.,  pg  301 , 
Wylie  1 992 


Intergra! 


Like  on  like  /  ,niecter  '3ce 
impinging  orifices 


Cross  section  through  an  etched  plate 


Oxidizer  distribution 
manifold  plate 


Four  sOcked  plates 
(with  exaggerated  thickness) 


Enlarged  detail  segment  of 
one  design  of  iike-oivlike  impinging 
orifices  and  feed  passage  geometry 


Separate  slides 


®  SSME  cutaway 
©  SSME  injector 

•  SSME  single  element  injector 

•  Platelette  injector  concept 
©  Henken  platelette 

©  Spashplate  platelette 


Overall  characteristics 


'  Pressures  (main  chamber): 

500  -  4000  psi 
1  Flow  rates: 

■  very  small  (satellite  thrusters)  to 
rates  that  can  drain  average 
swimming  pools  in  seconds 
(boosters). 

Combustion  chamber  volumes: 

■  0(1  in3)  (satellite  thrusters)  to 
0(1  ft3)  (boosters) 

Propellants: 

■  GOX,  LOX,  H2,  RP-1  (kerosine), 
NTO/MMH,  many  more. 

Single  element 
e—10  Ib/s/element 


\\  V.  V.;-.  I 

'  '  ; 

AFRL.  test  stand 
1A  flare  stack  test 


•  -V>  •  •  ?pj 
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.  .  .v."  ■  ‘  1 


Air  Puree  Rcseiirih  Lull  oratory 
P r i >j )  u Isi on  Direct oni (c 
Edwards  Research  Site 

Test  St;uid  I A  is  a  Ij^uiel  rocket  engine  test  stand  designed  Tor  static  firing  the  largest  rocket  engines 
ever  built,  with  thrusts  up  to  1.600,000  pounds.  Originally  built  in  1956  for  the  Atlas  Intercontinental 
Ballistic  Missile  Program,  the  stand  was  modified  into  a  rcx±et  engine  st:uid  in  I960  for  die  Apollo 
program.  Tile  first  stage  engine  of  the  Saturn  V  moon  rocket,  the  F-l ,  was  test  fired  on  1 A  for  most  of 
the  Apollo  years.  In  1995,  a  decision  was  made  to  modify  Test  Stand  1A  for  the  Evolved  Expendable 
Launch  Vehicle  (EELV)  program.  Liquid  hydrogen  systems  were  added  to  the  stand  to  provide  the 
capability  of  testing  cryogenically-fueled  engines.  The  instrumentation  and  control  systems  for  the 
stand  have  been  completely  upgraded  to  the  latest  state-of-the-art. 

The  IA  test  stand  is  tut  integral  part  of  the  1-120  Large  Rocket  Engine  Development  and  Test  Facility, 
where  rocket  engine  technologies  for  the  future  are  being  developed  for  the  U.S.  Air  Force.  Testing 
capabilities  have  been  improved,  allowing  both  commercial  and  Government-funded  programs  to 
coexist  within  the  test  facility,  maximizing  useof  the  facility  resources. 

Test  Stand  1A  Specifications: 

90,000  gal  Liquid  Hydrogen  Run  Tank  1.6Mlbf  Thrust  Measurement  System 

75,000  gal  Liquid  Oxygen  Run  Tank  60,000  HP  Facility  Turbine  Drive  System 

300  Ib/see  Hydrogen  Flare  Capability  Separate  Kerosene  Run  T;mk 

320  channel,  100K  sample/sec  Digital  Data  Acquisition  System 
Flame  deflector  with  1.000,00  gallons  cooling  water 


Main  Photo:  Test  Stand  1A  during  a  liquid  hydrogen  flow'  test.  Tire  large  flame,  which  peaked  at 
almost  700  ft  above  grade,  was  the  result  of  flowing  200  ib/sec  liquid  hydrogen  for  5  seconds  into  the 
main  flareslack.  The  1A  superstructure  is  100  feel  high.  The  structure  to  the  right  in  the  foreground  is 
Test  Stand  2A,  the  facilities’  component  development  stand. 

Left  Inset:  Test  Stand  1 A  during  the  installation  of  the  90,000  gallon  liquid  hydrogen  run  lank.  The 
900  ton  Manitowoc  crane  is  one  of  the  largest  mobile  cranes  in  the  world.  The  facility  is  positioned 
over  the  side  of  a  mountain,  allowing  the  rocket  exhaust  to  travel  150  ft  downward  before  impacting 
the  water-cooled  flame  deflector. 


Right  Inset:  Tills  1997  ribbon-cutting  ceremony  shows  the  completed  lest  stand  ;uid  gives  a  si 
perspective  of  the  facility. 


•  Injector  requirements 

i  ^  .... 


<»  Complete  combustion  in  the 
shortest  possible  length 
a  Main  injectors: 
performance  vs  weight 
tradeoffs 

0  Preburners/GG’s: 
downstream  component 
interactions,  eg,  turbine 
blades,  etc 
»  Acoustically  stable 
0  Chamber  modes 
b  Feed  system  coupling 
•  Chamber/wall  compatibility 
b  Heat  transfer/cooling 
■  Oxygen  blanching 
b  Lifetime 


•  Minimize  pressure  drop 
.*  Throttling 

•  Ignitable;  minimum  ignition 
transients 

•  Cost,  weight 

•  The  “ilities:” 

■  Reliability 
b  Maintainability 
b  Manufacturability 
b  Durability 
b  Operability 
b  PREDICTABILITY 


Major  trends  since  the  end  of  the 
cold  war 


m 


Major  trends  since  the 


end  of  the  cold  war 


•  Infusion  of  Russian  (formerly  Soviet)  technology 

•  Continued  rise  of  cost  as  a  major  consideration. 
Injector  impacts: 

h  Aversion  to  risk 

a  Interest  in  using  CFD  tools  combined  with  subscale 
data 

•  Gas/gas  injectors 

©  Understanding  of  injection  phenomena  at 
supercritical  pressures. 


Russian  experience:  injector  impacts 

•  fc- _ _ _ _ _ „ _ _ _ _ 

mfm 

Relative  Soviet  disadvantage  in  machining 

&  Fewer,  coarser  elements,  larger  thrust-per- 
element 

■  0(1 0,000- Ibf/element)  vs  0(1,000)  Ibf/element  in  the 
US. 

©Tendency  towards  coaxial  elements,  typically 
swirl  coaxial 

■  Counter-swirl  variants 

©Injectors  were  highly  performing  despite 
coarseness 


Russian  experience:  injector  impacts 


"21 

Use  of  coatings 

®  Injector  impact:  mitigates  against  streaking 
caused  by  coarse  elements 

Use  of  oxidizer-rich  preburners 
©  Partially  enabled  by  use  of  coatings 
©  Contrary  to  conventional  US  wisdom  at  the  time 

•  Contributed  to  US  interest  in  full  flow  staged 
combustion  cycles 


staged  combustion  cycle 

Advantages 

®  Eliminate  F/O  seal  in  L02 
pump 

■  Reduce  weight 

■  Eliminate  catastrophic  failure 
mode 

»  Reduce  Ox  turbine  input  temp 
•  Reduced  maintenance  and 
failures 

*  No  cryogenics  downstream  of 
preburner;  “standard”  piping 

Injector  impact 

*  Gas/gas  main  injectors 

■  No  previous  experience  at 
main  injector  scales 


Injector-related  developments  of  the  1990’s 


Kistler  K-1 
Long  March 
Rotary  Rocket 
Scorpion 
Shuttle 
Titan  IV 
X-33 
X-34 


2002  RD-0180 

2002  RL-tOA-2 
1S36  Vulcain 
1 996  Aeslus 
dev  VINCI 
dev  Beal 

RL-10A-5 

RS-68 

RL-10-B-2 

IE-7 

LE-5E 

NK33.  Nk43 

YF-40 

RR 

Scorpion 
SSME  block  II 
LE-87-1 1,  LE-91-1 


Faslrac 
dev  dev 

RL- 50/60 

Integrated  Powerhcad  Demo 
TRW  600  Klb!  pintle 


Stacie  Manut. 

main  NPO  Encrgomash 

upper  Prall 

main  SEP 

upper  DaimlerChrysler 

upper 

main  Beal 
main  Pratl 
main  Rocketdyne 
upper  Prall 


main  RR 
main  Microcosm 
main  Rocketdyne 
Aerojet 

main  Rockeldyne 
main  NASA 
dev  dev 
upper  Pralt 

boost  Rockeldyne/Aorojot 
boost  TRW 


1 999  Lox/kerosine 

1995  Lox/H2 

1996  Lox/H2 
1996  N204/MMH 
dov  Lox/H2 

dov  H202/kerosine 
Lox/H2 
Lox/H2 
Lox/H2 
Lox/H2 
Lox/H2 
Lox/kerosine 
N204/MMH 
H204/korosine 
Lox/kerosine 
Lox/H2 

N204/Aoro2inc  - 
dov.  Lox/H2 

Lox/kerosine 
dev  dov 

Lox/H2 

Lox/H2 

Lox/H2 


Rockot  Propulsion  Technology  program 


Twochamb.  ver.  RD170  (1970’s) 

Derivative  of  1 95Q’s  engine 

Europe 

Europe 

Europe 

Privately  funded 

EELV 

EELV 

Japan 

Japan 

Russian  engine  on  US  vehicle 
China 

Privately  fundod 
Privately  funded 
upgrade 


NASA  Space  Launch  Initiative 


Injector-related  developments  of  the  1990’s  + 


General  Trends 

Risk-adverse  tendency  to  stick  with  proven 
injector  types,  with  some  exceptions: 

■  Integrated  Powerhead  Demo  preburner  and  main 
injectors 

*  Gas/gas  main  injector  development  (discussed  later) 

*  Hydrogen  peroxide  injectors 

*  Nearly  all  proprietary -not  discussed  here 

■  Notable  mention  given  to  pintle  injector  development 

&  Although  the  basic  type  was  not  changed,  there  was  a 
respectably  large  extrapolation  beyond  existing  experience  (to 
a  600  Klbf  H202  test). 


Recent  developments  in  liquid  rocket 
injector  design  tools 

Outline 

©  Overview  of  atomization 
©  Recent  developments  in  modeling 
®  Recent  developments  in  experimental  methods 
h  Cold  flow  characterization  of  rocket  injectors  at  AFRL 


Overview  of  atomization 


Atomization  steps 


Lass 


Primary  atomization 


Impinging  element  spray;  plane  of 
jets  perpendicular  to  slide 


Helical  structures 
in  pre-impingement 
jet 


*  Vaporization/combustion 


Breakup  of  jets  and  sheets  into  Gasification  so  that  fuel  and 


long  irregularly  shaped  ligaments 
Secondary  breakup 
Breakup  of  ligaments  into  droplets 
Breakup  of  drops  into  smaller  drops 


oxidant  can  mix  at  the 
molecular  level  and  burn 


Breakup  mechanisms 


©  Shearina  Stretch-Thin-Shear  illustrated  for  a 

y  typical  prefilming  atomizer  (non- 

©  Liquid  phase  turbulence  rocket) 

©  Surface  tension  LIQUID 

,jL  '  PREFILMIh’G 

An  efficient  atomization  0  . 

mechanism:  F^-l 

Stretch-Thin-Shear  r 

_  C  y/y 

(stretch  liquid  into  the  ^ 

thinnest  possible  sheets,  Ymrrm / / rrm / b pfa 

then  shear)  /  /ATOMtItNO  E3GE 


_  ATOMIZING  EDGE 


Rizk,  N.K.,  and  Lefebvre,  A.H.,  ‘The  Influence  of  Liquid 
Film  Thickness  on  Airblast  Alomizalion,”  Journal  of 
Engineering  for  Power,  vol.  102,  pp.  706-710,  July  1960. 


Overview:  drop  size  distributions 


»  Nukiyama  and  Tanasawa  (four  parameter) 
dN 

_ zz  ciDP  exp [—(bD)°]  experimental  distribution 

clD  (form  unspecified) 

where  a,b,p,q  are  parameters 


»  Rosin-Rammler  (two  parameter) 

1  —  Q  —  eXp[— (Z)  /  X y1  ]  Particle  Diameter  (//m) 

where  Q  is  the  fraction  of  the  total  volume  contained  in 
drops  of  diameter  less  than  D,  X  is  the  diameter  such  that 
63.2%  of  the  total  liquid  volume  is  in  drops  of  smaller 
diameter,  and  q  is  a  parameter. 

©  Several  others 


Overview:  representative  drop  diameters 


Mean  diameters 


V(a-b) 


a  -V  b  Name  of  mean 

a  b  (order)  Symbol  diameter 


Expression  Application 


1  0  1 

2  0  2 

3  0  3 

2  1  3 

3  1  4 

3  2  5 

4  3  7 


Dl0 
D » 
Ao 

A , 

Z>J, 

A: 

A> 


Length 
Surface  area 


Volume 


Surface  area-length 
Volume-length 

Sauter  (SMD) 

De  Brouckere  or 
Herd  an 


INfr 


2A£? 

2NiDl 

/XNfDlV'2 

XSN'dJ 

IV'D* 

2N,D] 


Comparisons 

Surface  area 
controlling 

Volume  controlling, 
e.g.,  hydrology 

Absorption 

Evaporation, 
molecular  diffusion 

Mass  transfer, 
reaction 


Combustion 

equilibrium 


Overview:  representative  drop  diameters 


Median  diameters:  Dr 


where  xx%  of  the  total  liquid  volume  is  in  drops  of 
smaller  diameter 

B  A).5  is  known  as  the  Mass  Mean  Diameter,  MMD 

-NOTE:  when  the  shape  of  the  distribution  is 
known,  it  is  only  necessary  to  report  one 
diameter,  plus  ratios  of  diameters  equal  to  the 
total  number  of  free  parameters  in  the 
distribution. 


Rocket  mean  diameters  from  the  1960’s 


®  Tend  to  be  reported  as 

D 43 > 

D  -5>'D'4 

43 "  5>A3 

which  is  a  volume- 
weighted  mean  diameter 

®  Reason:  old  data  used 
the  molten  wax  technique 

a  Sieving  process  gives  a 
volume  weighted  mean. 


BUI 


mmmm 


Wax  sieve  cuts 


Sample  advancement:  secondary  breakup 


Secondary 

breakup 

regimes 


n-HEPTAMC  o  9LYCER0L 

□  GLYCEROL 
OLY CEROL  2  I X  a  GLYCEROL 

OLY  CEROL  63%  -  UERCUwr 

GLYCEROL  7fl%  O’  eoOFLUlO 

8LYCER0L  14%  » 


GLYCEROL  62%  »  HWSW  ,f  eL  UM3J  6 

GLYCEROL  87%  4  WINZE  (1053)  * 

GLYCEROL  »afl%  r  LANE  ( 1901)  4 

MERCURY  ♦  LOf>AREV(»T3[  « 

COOFLUlO  ►  KRZecZKOWBKIll9B5) - 

PRESENT  - 


Pi^Po  .  8*0-12000 


We  =  p  d0u o  /  a 


Oh  =  n,  /(p,d0o)u 


Chou,  Hsiang  and  Faeth, 
Int.  J.  Multiphase  Flow,  Vol. 
23, 1997 


10%  <  DEFORMATION  <  20% 
D%  <  DEFORMATION  <  IO%  , 


Ay 


•-/// 


id*  to0  \6 
OHNESORGE  NUMBER 


,(r  Advancements  in  secondary  breakup 

ifi“  ~  ™  . . . .  •  -4 . 

©  Regimes  mapping  when  different  kinds  of 
breakup  occur,  but  that  was  all 

©  Chou  and  Faeth  correlated  data  for  the  bag 
breakup  regime  which  gives  a  complete  picture: 

■  How  long  it  takes  to  breakup 

“  Size  of  the  parent  droplet  when  breakup  stops 

■  Size  distribution  of  all  the  daughter  droplets 

©  NOTE:  The  results  of  this  sample  advancement 
pertain  to  liquid-to-gas  density  ratios  of  about 
500,  much  larger  than  in  most  rockets 

m  Many  recent  advancements  in  atomization  don’t  pertain 
to  rocket  conditions 

Chou  and  Faeth,  Int.  J.  Multiphase  Flow,  Vol.  24,  1998 


Recent  developments  in  modeling 


^  Computer  resources 

|ll  ■  - - - - 

©  Increases  in  processor  speed 
a  Currently  0(2.2GHz) 
h  <1 00MHz  1 0  years  ago. 

©  Massively  Parallel  Computing 
■  Multiple  processor  computers 
m  Clusters 
e  Reduced  cost 

©  Increased  computer  power  provides  ability  to 
compute  flows  with  more  physics,  and  with 
greater  resolution  and  speed. 


Advances  in  communication  between  codes 


®  Typically  the  design  of  components  such  as  thrust 
chambers  involves  the  use  of  more  than  one  code  to 
perform  the  required  calculations 
»  Feed  system  codes  ■  Film  cooling  codes 

■  Injector  codes  ■  Nozzle  codes 

b  Combustion  chamber  cooling  codes  B  Instability  codes 

®  Use  of  different  codes  requires  transferring  information 
between  them,  which  can  be  extremely  tedious 

®  Executive  codes  are  being  developed  which  automate  the 
hand-offs  between  codes,  saving  significant  time 

■  Boeing  Thrust  Chamber  Analysis  Toolkit  (TCAT) 


Advances  in  gas-side  CFD 


•  More  robust  numerical  algorithms 

b  Preconditioning 

Choi,  Y.H.  &  Merkle,  C.L.,  “The  Application  of  Preconditioning  in  Viscous 
Flows,”  J.  Comp.  Phy.,  72,  1987. 

b  Discretization  schemes 

*•  Multigrid 

Mavriplis,  D.J.,  “Multigrid  Techniques  for  Unstructured  Meshes,”  ICASE  Report 
No.  95-27,  April  1995. 

«  Conservation  Element/Solution  Element  (Wang  et  al.) 

Wang,  X.-Y.,  Chow,  C.-Y.,  &  Chang,  S.-C.,  “Application  of  the  Space-Time 
Conservation  Element  and  Solution  Element  Method  to  Two-Dimensional 
Advection-Diffusion  Problems,"  NASA-TM-1 06946,  1995. 

vAdvection  Upstream  Splitting  Method  (AUSM) 

Liou,  M.-S.,  &  Steffen,  C.J.,  “A  New  Flux  Splitting  Scheme,”  J.  Comp .  Phy., 
107,  1993. 


^  Advances  in  liquid-side  GFD 


©  Surface  Tracking  Algorithms 

a  Volume  of  fluid 

Bussmann,  M.,  Mostaghimi,  J.(  &  Chandra,  S.,  Phys.  Fluids,  11,  1999. 

■  Free  surface  tracking 

Helenbrook,  B.T.,  Comp.  Meth.  Appl.  Mech.  Eng.,  191,  2001. 


CFD  for  rocket  applications 


•  Improved  computer  power  allows  more  physics  to  be  included  with 
greater  resolution,  but  the  physical  models  continue  to  remain 
mostly  not  validated  for  rocket  conditions 

■  High  pressure  and  temperature 

b  Supercritical  fluids 

h  Dozens  of  species  and  scores  of  reactions 

•  Computer  resources  are  still  not  adequate  for  simulating  complex 
injector  phenomena 

■  Atomization  is  sensitive  to  upstream  manifolding;  upstream 
manifolding  must  also  be  modeled 

a  Modeling  atomization  requires  resolving  turbulent  length  scales 
inside  tiny  orifices 

■  Computing  rocket  design  parameters  such  as  wall  heat  fluxes 
requires  resolution  of  the  entire  combustion  chamber 


rocket  applications 


Consequence 

©  Injector  atomization  performance  currently  tends 
to  be  an  experimental  input  into,  not  a 
computation  resulting  from,  most  liquid  rocket 
combustion  design  codes. 


Recent  developments  in 
experimental  methods 


difficulty  /  cost 


The  problem 


•  Most  of  engine  development  cost  has  historically 
been  spent  on  trial-and-error  fixes  of  problems 
developed  after  full  scale  design  is  complete. 

©  Problems  not  discovered  until  full  scale  testing 
tend  to  be  extremely  expensive  to  fix,  and  have 
historically  required  sacrificing  original  engine 
performance  and/or  lifetime  goals. 


Hierarchy  of  injector  experiments 


General  trend:  an  increase  in 
subscale  efforts 


Injector  design 


Subscale 
Multi-element 
hot  fire 


Single  element 
'high  pressure; 
cold  flow 

Single  element;  . 

atmospheric 
.  '  cold  flow. 

Single  element 
. ;  hot  fire 
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degree  of  simulation 
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AFRL  cold  flow  facility 


What  good  are  subscale  experiments? 


The  performance  of 
these  two  injectors  is 
projected  to  be  essen¬ 
tially  equivalent 

*  Choose  the  one  that’s 
cheaper  or  that  has  some 
other  advantage 
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Thoughts  on  scaling 


©  Trend  analysis:  the  expectation  that  trends,  if  not 
magnitudes,  are  often  preserved  between  scales. 

■  Supported  by  experience  in  rockets  and  elsewhere  that 
injectors  which  do  relatively  better  in  cold  flow  tend  to  do 
relatively  better  in  hot  fire  -  eliminate  obvious  bad  designs. 

©  Bracketing  and  limiting:  case-dependent  projections 
of  the  directions  in  which  performance  at  one  scale 
will  deviate  from  performance  at  another  scale. 

■  Ex:  an  injector  that  has  a  c*  efficiency  of  95%  as  a  single 
element  may  often  be  expected  to  have  a  c*  efficiency  of 
95%  or  better  as  multi-element. 

h  Ex:  if  you’re  burning  something  up  at  low  pressures, 
chances  are  you’re  going  to  burn  it  up  at  high  pressures. 


.  Thoughts  on  scaling 

I  . . . ,. . . . . 

e  Code  validation:  if  validated  at  one  scale,  codes  may 
often  have  the  necessary  correct  physics  to 
extrapolate  reliably  to  another  scale. 

©  Physical  understanding:  If  mysterious  things  are 
happening,  any  scale  which  provides  understanding 
is  usually  useful. 


L  Diagnostics 

St- 


At  almost  any  scale,  optical  diagnostics  remain 
notoriously  difficult. 
b  Often  impossible  at  too  large  a  scale 
Prior  to  the  1990’s: 

■  Cold  flow  injector  data  was  limited  to  about  400  psi  or 
less 

■  Hot  fire  data  was  limited  largely  to  simple  visualizations 
At  present: 

bAFRL  routinely  characterizes  injectors  in  cold  flow  at 
pressures  up  to  2000  psi. 

m  Many  examples  of  advanced  diagnostics  applied  in  hot 
fire  exist 


Diagnostics 


The  development  of  effective  subscale 
methodologies  has  been  as  much  about  the 
development  of  diagnostics  that  work  as  it 
has  been  about  injectors 


»  Most  existing  diagnostics  have  been  developed  for  low 
pressure  applications  and  have  significant  difficulties  at 
high  pressures. 

*  Much  if  not  most  of  the  effort  in  the  1990’s  was  spent 
developing  diagnostics  for  only  one  kind  of  injector  - 
coaxial  injectors 


Characterizing  rocket  injectors  at 
AFRL 


,  Uses  of  single  element  rocket  injector 
!j t  measurements 

iiKti.  mi* 

•  Understand  physics 

©  Provide  input  to  existing  codes  that  require 
injector  performance  as  an  input 
s  Need  accurate  data 
®  Validate  codes 
e  Need  accurate  data 

©  Evaluate  particular  injector  designs 

■  Not  providing  an  answer  is  the  same  as  letting  the 
designer  go  with  his/her  existing  superstitions 

a  Need  to  consider  the  data  together  with  all  its 
uncertainties  and  still  give  the  best  possible  guidance 


AFRL  full  scale  single  element  cold  flow 
injector  characterization  facility 


@  Atmospheric  tests  are  used  to: 

■  Check  for  manufacturing  defects 
h  Check  out  diagnostics 
®  Elevated  chamber  pressures  are 
required  to: 

h  Prevent  cavitation  when  at  realistic 
pressure  drops 

a  Match  more  scaling  parameters 
©  AFRL  facility: 

e  Chamber  pressures  to  2000  psi 
using  water  as  a  simulant. 
h  Mechanical  patternation  at  full 
pressure. 

*  Full  suite  of  optical  diagnostics 


ram 


Relative  Mass  Fraction 


Pressure  capability  matches  more  scaling 
parameters 

SsfclsIfcU  - 


Atmosph eri c  Pressure 
Geometry  Scaling=1.0 
Velocity  Ratio  O/F=.084 
Ma=0.25 


Pc=43j>  mi  g 
Geometry  Scaling=1.0 
Velocity  Ratio  O/F=.084 
Ma=0.25 

Momentum  Ratio  O/F=.079 
Density  Ratio  0/F=29.2 


SS'ME  fuel  preburner  coax  injector  tests  demonstrate 
reduced  effectiveness  of  swirl  on  macro-scale  mixing  at 
pressure. 


Shear 


First  step:  always  take  a  picture 


Macro  (shadow¬ 
graph,  Schlieren, 
laser  sheet) 

*  Gross  features 

■  Relative  fineness 
h  Plan  strategy 

■  Videos  usually 
sufficient 

Micro  (sometimes) 
b  Drop/blob  sizes 


Sprays  are  almost 
always  very  dense 
«  Light  extinction 
often  reaches 
100% 


:  Second  step:  mass  distributions 

m - - - . — . 


«  AFRL  uses  a  traversable  linear 
mechanical  patternator  to  measure  2D 
mass  flux  distributions  at  pressure. 

*  Measuring  collection  efficiencies  a  must, 
m  Liquid/liquid  injectors:  ?7C0|,  =  80-90% 
b  Gas/liquid  injectors  with  large  gas 
momentum:  ?jcoi,  =  30-50% 
b  Current  best  practice  is  to  correct 
patternation  results  by  the  collection 
efficiency 

<*  In  most  cases  we  project  that  intrusive 
errors  cause  measured  profiles  to  be 
smoother  than  they  actually  are. 


Optical  patternation  methods 


Illumination  of  the  spray  with  laser 
sheets 

«  Of  qualitative  value,  but  prone  to 
quantitative  errors  in  dense  sprays 


*  Extinction  of  the  laser  light 
■  Attenuation  of  the  signal  light 

*  Secondary  scattering  from  particles  outside  the 
plane 

©  Significant  advances  in  developing 

quantitative  corrections  were  made  in  the 
1990’s,  but  these  tend  to  remain 

insufficient  for  rocket  applications 

Errors  caused  by  light  ekiriction 
(Illumination  from  left) 


4?  Thirc*  steP:  drop  sizing  and  velocimetry 

is*  &&  - 

Categories  of  optical  techniques 

e  Imaging  (digital,  film) 

■  Image  processors  can  perform  focus  discrimination  by  measuring 
edge  gradients,  and  automatically  measure  sizes. 

a  Depth-of-field  depends  on  size;  large  drops  preferentially  sampled 
?  Corrections  can  be  made 

■  Main  advantage  is  ability  to  measure  non-spherical  droplets 

■  Common  in  rocket  sprays 

■  Main  disadvantage  is  slow  sampling  rates 

*  Few  significant  advances  in  the  1990’s  except  for  faster  computer 
speeds  and  increased  resolution  of  CCD  cameras. 

»  Phase  Doppler  interferometry  (covered  next) 

«>  Forward  scattering  methods  (covered  next) 


Phase  Doppler  interferometry 


Droplet  Size  and  Velocity 

•  Spatially  and  temporally 
resolved. 

■  Single  droplet  counting 
technique. 

■  Extension  of  the  laser 
Doppler  velocimetry 
technique. 

u  Size  range  from  0.5  to 
1000  pm  (droplets  must 
be  spherical). 

■  Up  to  three  velocity 
components,  droplet 
size  and  mass  flux. 


Transmitter  Lens 
Beamsplitter  j  j 

Laser 


I  ^——4  Receiver  Lens 


jfc'l  Frequency  and  Phaser 
f,  \  Sigfial  Processor  . 


8acha!o.  Alomiialion  and  Sprays,  Vol.  10,  2000 


Phase  Doppler  interferometry 

Simultaneous  measurement  of  droplet  size  and  velocity 

*  Phase  shift  proportional  to  droplet  diameter. 

23  3  detectors  yield  2  phase  measurements. 

Redundancy  check 
Prevents  ‘‘wraparound” 

-  Sphericity  check 

■  160  MHz,  1-bit  sampling,  Fourier  transform  “real-time”  processing. 


BachaloandHousor.  AIAA-84-1 199,  t$aj  Fijurf  7.  Phue/Dopple»  initrumtnt  response  curte 


Phase  Doppler  interferometry 


Tracking  “gas”  and  liquid  phases 

*  For  gas/liquid  injectors,  accounting  for  both  phases  is  important 

*  With  PDI,  the  smallest  drop  sizes  can  be  assumed  to  mainly  follow 
the  gas.  Therefore,  filtering  for  the  smallest  drops  traces  the  gas  flow 

■  Allows  estimation  of  the  mixture  ratio 

Non-spherical  particles 

•  Size  measurements  don’t  work  for  non-spherical  drops,  and  non- 
spherical  drops  are  not  always  rejected  by  the  signal  processor 

“  Photographic  inspection  of  the  measurement  location  should  be 
performed  to  ensure  that  most  of  the  droplets  are  spherical 

•  Valid  velocity  measurements  are  still  obtained  by  operating  the 
instrument  as  an  LDV. 


Phase  Doppler  Interferometry 


Limitations  in  dense  sprays 

a  Designed  for  “dilute”  sprays. 

a  Requires  single  droplet  in  the  probe  volume. 

a  Beam  waist  usually  much  larger  than  maximum  droplet  size 
(uniform  illumination). 

«  Dense  sprays  result  in  multiple  droplet  occurrences,  beam 
attenuation  and  multiple  scattering  of  signals. 

0  1990's  development  -  Reduce  probe  size  to  be  much  smaller  than  the 
largest  droplet  to  be  measured. 


Technique  works  because 
only  a  fraction  of  the 
droplet  surface  is  actually 
imaged  into  the  receiver 


Area  of  Particle 


w  Phase  Doppler  interferometry 


Small  Probe  Volume  Technique 

•  Solving  large  probe  volume  problem  creates  another: 
Trajectory  dependent  scattering  errors. 

■  Uneven  illumination  can  result  in  significant  amounts  of 
reflected  light  reaching  the  receiver. 

«*  Geometric  optics  modeling  indicates  that  these  errors  can 
be  rejected  with  intensity  validation. 
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Small  probe  volume 

technique  is  ineffective 
when  light  extinction 

reaches  100%. 

Preferred  approach  is  to  use 
beam  guides  /  flow  splitters. 
*  Not  always  possible. 
h  Flow  effects  assessed  by 
noting  trends  as 
separation  distance  is 
varied. 


Splitter  Plate  Separation  (mm)  Splitter  Plate  Separation  (mm)  Splitter  Plate  Separation  (rrm) 


,  „  Phase  Doppler  interferometry 

IIS . . ~ . . . 

Small  probe  volume  technique 

®  More  complex  than  conventional  PDI. 

■  Requires  simultaneous  measurement  of  phase  and  peak  scattered 
light  intensity. 

®  Small  beam  waist  (~  60  pm)  results  in  short  transit  times  requiring 
high-speed  data  sampling  (160  MHz). 

■  Also  requires  a  new  method  of  probe  volume  correction  (PVC)  for 
accurate  droplet  sizing  and  mass  flux  measurements. 

•  Combined  with  a  flow-splitter,  the  small  probe  volume 
technique  can  greatly  improve  accuracy  of  mass  flux 
measurements  in  sprays  with  droplet  number  densities 
approaching  105  droplets/cm3. 


Flux(cc/s/cm  2) 


Phase  Doppler  interferometry 

Effect  of  small  probe  volumes  and  flow 
splitters  on  measured  mass  fluxes 


Conventional  Probe  /  No  Splitter 
Conventional  Probe /With  Splitter 
-re—  Small  Probe/ With  Splitter 
Patternator  * 


SSMB  Prchnrner  Injector 
Pc=2.9  MPa 
f-T20=32.7  gm/s 
N2=35.3  gm/s 
Z=50  mm 


Radial  Position  (mm) 


Forward  scattering  instruments 


«  Ensemble  forward  scattering 

a  Yields  droplet  size  distribution. 

■  Transmissions  down  to  2%  with  multiple  scattering 
correction. 

■  Size  range  from  0.1  to  2000  (xm  (forgiving  of  slightly  non- 
spherical  droplets). 

h  Droplet  distribution  calculated  from  inversion  of  light 
scattering  distribution  on  a  31  ring  detector. 


small  drop 


Ring 

Detector 


large  drop 


0=  scattering 
angle 


Beam  Power 
Detector 


’m  Forward  scattering  instruments 

Ring  Detector  Approach 

•  Detector  measures  the  diffraction  pattern  imaged  by  the 
Fourier  transform  lens. 

■  Small  droplets  have  large  scattering  angles, 
a  Large  droplets  have  small  scattering  angles. 


Beam  Power 
Detector  N 


Sea  tiering 
Detector 


a,  .  , 
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Beam  Focus 
&  Pin  Hole 


Mie  Scattering  Calculations 
n  =  1 .33  -  O.Oi 


0.00  0.05  0.10  0.15  0.20  0,25  0.30 

Scattering  Angle  (0) 


Forward  scattering  instruments 

Scattering  Theory 

•  Earlier  instruments  used  Fraunhofer  diffraction 

theory. 

a  Limited  to  droplets  much  larger  than  the  laser 
wavelength. 

m  Fixed  scattering  cross-section  (C  =  2  *  A  .  ) 

a  '  cross-section/* 

m  Does  not  account  for  anomolous  diffraction  (reflection  + 
refraction). 

®  Newer  instruments  (1990’s+)  use  Lorenz-Mie 

theory. 

a  Mostly  due  to  faster  computers. 

a  Eliminates  limitations  imposed  by  Fraunhofer  diffraction 
theory. 


w  Forward  scattering  instruments 

Recent  Advances 

a  Temporally  resolved  data  (faster  computers). 

®  Use  of  Mie  scattering  as  opposed  to  Fraunhofer 
diffraction  theory. 

®  Better  multiple  scattering  corrections  at  large  beam 
obscurations  (model  independent  analysis). 

Limitations 

»  Moderately  dense  sprays  (Transmission  >  2%). 

*  Can  be  sensitive  to  refractive  index  gradients  (high- 
pressure  &  turbulent  beam  steering). 

*  Depends  on  drop  sizes 

*  Poor  spatial  resolution  (line-of-sight,  10  mm  beam  diam.) 


Forward  scattering  instruments 


Dense  spray  corrections  are  used  with  forward 
scattering  techniques 


Regime  often  of 
interest 


•  Stirred  cells  of  polystyrene  beads 
currently  being  used  to  assess 
potential  magnitude  of  errors. 

«  Approaches  to  index  of  refraction 
gradients  masquerading  as  large 
droplets  also  being  investigated. 

Glass  Cell  with  Known 
Distribution  of 


Magnetic  Stirrer 


Forward  scattering  instruments 


Beam  wandering  effects  at  high  pressures 

©  Index  of  refraction  gradients  cause  beam  to 
wander  off  the  center  hole,  hitting  inner  rings  and 
appearing  to  be  large  droplets 


II 

If 

No  saturation  of  gas  with  water  vapor 
Temperature  control  to  +/-  7.0  C. 


Saturation  of  gas  with  water  vapor 
Temperature  control  to  +/-  0.7  C. 


Forward  scattering  vs.  PDI 


Do  forward  scattering  instruments  and  PDI 
give  the  same  answer? 

©  Generally,  no. 

n  PDI  is  a  point  measurement  (where  the  beams  cross), 
while  forward  scattering  is  integrated  along  the  length 
of  the  droplets  in  the  beam. 

a  PDI  is  flux-based,  while  forward  scattering  is  volume- 
based 

•s-As  a  point  measurement  recording  droplets  going  by  a  fixed 
location  in  space,  PDI  is  biased  towards  size  classes  that  are 
the  fastest  moving 

Forward  scattering  is  biased  towards  size  classes  that  have  the 
greatest  number  density  in  the  beam  at  any  given  time. 


^  Hot  fire  measurements 

•  Assuming  windows  are  available  in  the 
combustion  chamber  (a  non-trivial  accomp¬ 
lishment),  the  same  three  steps  generally  apply, 
except: 

a  Many  diagnostics  become  significantly  harder  in  hot 
fire;  consequently,  the  range  of  choices  becomes  more 
limited. 

©  One  silver  lining  is  that  hot  fire  tends  to  quickly 
burn  off  the  fine  particles, 
a  Optical  obscuration  is  sometimes  improved 


If  Hot  fire  measurements 

®  Measurements  that  have  been  relatively 
“straightforward”  (i.e.,  difficult) 

■  Visible  emission 
a  Shadowgraphy 
a  Schlieren 

a  Emission  spectroscopy 
a  Line-of-sight  absorption 
«  Disadvantages 

■  All  these  are  line-of-sight 
'■t  Beam  steering  effects 


w  Hot  tire  measurements 

-j  . . . . 

Measurements  that  have  been  achievable 
with  significant  difficulty 

•  Velocimetry  (laser  Doppler  velocimetry  (LDV), 
particle  imaging  velocimetry  (PI V)) 

•  Phase  Doppler  interferometry 
©  Mie  scattering 

■  Location  of  liquids 


<ajj.  Hot  fire  measurements 

Measurements  that  have  been  achievable 
with  great  difficulty 
®  Spontaneous  Raman  spectroscopy 
■  Species  concentrations 

©  Laser  induced  fluorescence  (OH,  CH,  etc.) 

“  Limited  quantitative  concentrations,  but  .  give 
reasonable  locations  of  the  flame  front. 

©  Coherent  anti-Stokes  Raman  Spectroscopy 
(CARS) 

m  Species,  temperatures 


m 


Major  trend 


Spontaneous  Raman  spectroscopy 


•  A  weak  effect  where  inelastic  scattering  causes  a 
shift  in  wavelength  by  an  amount  dependent  on  the 
molecule. 


ta  Used  primarily  for  species  concentration  measurements 
©  Signal  strength  proportional  to  density,  inversely 
proportional  to  the  fourth  power  of  the  wavelength 
h  Favors  high  pressures  and  blue  and  UV  wavelengths 
h  Works  best  in  single-phase  flows;  multiphase  effects 
cause  interference. 


®  Most  of  the  major  labs  have  developed  this 
diagnostic. 


Case  study:  gas/gas  injector 
development 


¥m, 
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Gas/gas  injectors 


e  Two  reasons  to  study  gas/gas  injectors 

h  Motivated  by  US  interest  in  a  full  flow  staged  combustion  cycle 
■'  Required  a  gas/gas  main  injector  to  be  developed  at  a  size  for  which 
the  US  had  no  experience 
m  Logical  first  step  for  research 

•  Tendency  to  model  or  measure  gas/gas  first  before  introducing  the 
potential  complications  of  multiphase  flows 

*  The  Gas/Gas  Injector  Technology  program  (GGIT)  was 
initiated  in  the  mid  90’s  by  a  team  composed  of 
government,  industry,  and  university, 
o  The  effort  attempted  to  combine  data  obtained  at  different 
scales  with  CFD  to  develop  a  design  for  the  injector 

*  Tucker,  et.  al.,  paper  AIAA  97-3350. 


GGIT  objectives 


©  Provide  Rocketdyne  with  data  on  candidate 
gas/gas  injector  designs  for  its  planned  FFSC 
engine  (RS-2100) 

&  Develop  a  national  gas/gas  database 
©Validate  CFD  codes  for  future  gas/gas  injector 
design 


The  team  evaluated  designs  which  were 
proprietary  to  Rocketdyne  and  other  designs 
which  were  in  the  public  domain 


GGIT  team  responsibilies 


®  Rocketdyne  Propulsion  and  Power  (RPP) 

*  Provided  overall  technical  leadership 

n  Designed  and  fabricated  proprietary  injector  prototypes 
b  Internally  funded  its  own  participation 

®  Penn  State  (PSU) 

*  Evaluate  propellant  mixing  and  combustion  of  a  single  injector 
elements  (uni-element  testing) 

•  NASA  Lewis  (LeRC) 

■  Evaluate  propellant  mixing  and  combustion  of  a  multiple  element 
injectors  (multi-element  testing) 

•  NASA  Marshall  (MSFC) 

■  Managed  overall  funding  and  coordination 
b  Performed  CFD  modeling. 


GGIT  design  logic 


Figure  3.  Gas-Gas  Injector  Technology  Project  Logic  Diagram, 


ihe  Optically  Accessible  Chamber  at  PSU. 


©  The  program  was  discontinued  before  completion 
because  the  FFSC  concept  was  not 
downselected  for  RLV 
h  X-33  /  Venture  Star  was 

©  Several  proprietary  and  public  injectors  were  uni¬ 
element  tested. 

•  Only  the  proprietary  injectors  were  fired  as  multi¬ 
elements 
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GGIT  impact 


RPP  used  lessons  learned  in  GGIT  to  demonstrate  the 
combined  firing  of  an  ox  rich  preburner  and  a  gas/gas  main 
injector 

m  Ref:  Farhangi,  et.  al.,  paper  AIAA  99-2757. 

RPP  designed  the  main  injector  for  the  Air  Force’s  Integrated 
Powerhead  Demonstrator  engine  partially  using  GGIT  results. 
■  This  engine  is  still  under  active  development,  and  has  not  yet  been 
fired. 

A  data  set  survives  for  a  gas/gas  coaxial  element  that  has 
been  used  as  a  benchmark  case  for  model  development. 
m  Foust,  et.  al.,  paper  AIAA  96-0646. 


Gas/gas  CFD  modeling 


©  The  Penn  State  coaxial  gas/gas  data  has  been 
modeled  by  several  groups 
a  Merkle  code,  Penn  State 
b  AS3D  code,  DLR,  Germany 
h  FDNS  code,  MSFC 
b  CFD++  code,  AFRL 

»  Schley,  et.  al.,  paper  AIAA-97-3350'  , 

:  Archambault,  et.  al.,  papers  AIAA-2002-1088  "and  AIAA-2002- 
3594" 

©  All  exhibited  comparable,  although  different, 
agreement  with  the  data. 
h  Only  the  AFRL  results  will  be  presented  here 


Penn  State  configuration 


OH-radical  imaging 

Velocity  &  species  field  (Raman)  measurements 


Computational  configuration 


Solution  exhibited  extremely  slow  convergence  with  large 
residuals 

a  Due  to  the  attempt  to  force  a  steady  solution  on  an  inherently 
unsteady  flow 

®  Subsonic  injector  flow  meant  that  pressure  waves  inside 
the  injector  must  be  modeled 
®  Many  tricks  required  to  obtain  solution 
®  Objective:  compare  steady  solution  with  a  time-averaged 
unsteady  solution 


l-l,  Annulus 


Realistic  manifold  geometry 
made  no  difference 


Kadis!  Dist.  (in) 


/  Axisymmetric  results 


m. 


Time-Averaged  (last  6000  iterations) 
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Steady  (residual  ~10  2) 


Steady  solution  obtained  in 
-6000  iterations.  1 8,000 
time-accurate  iterations 
performed  after  steady  state 
reached. 

Independent  of  initial 
conditions  at  12000 

iterations. 

OH  concentration  contours 
demonstrate  differences 
between  steady  and 

transient  solutions. 

Steady  solution  does  not 
capture  the  effects  of  the 
flapping  flame. 


Comparisons  with  data 

H2  Mole  Fraction 


•  Gaseous  nitrogen  curtain  purge  in  experiment  not  modeled.  This 
accounts  for  decrease  in  H2  experimental  mole  fraction  near  walls  in 
the  data.  Possibly  also  accounts  for  computed  profiles  broadening 
more  rapidly  than  experimental  data 

»  Steeper  computational  profiles  may  be  due  to  unsteadiness  in  the 
experimental  shear  layer 


Comparisons  with  data 


02  Mole  Fraction 


Of  MoU  F»adion(X  =  2  ini 


Oxygen  does  not  radially  diffuse  as  does  the  hydrogen 
with  downstream  distance. 

02  and  H2  mole  fraction  calculations  agree  with  data 
about  as  well  as  other  reported  results. 


Comparisons  with  data 


Mean  Axial  Velocity  (m/s) 


©  Time-averaged  peak  velocities  better  match  data  than 
steady  velocities  near  the  injector.  Full  geometry  of 
injector  manifold  was  modeled  to  see  if  better  agreement 
could  be  obtained,  but  no  significant  improvement  was 
observed. 


Comparisons  with  data 


RMS  Axial  Velocity  (m/s) 


*  Steady  and  time-averaged  results  compare  about  equally 
well  with  experiment. 

®  Differences  in  inlet  turbulence  boundary  conditions  may 
likely  have  a  significant  effect  on  RMS  values. 


Effect  of  window  purge  (axisym.) 


Hj  Molo  Fraction  (X  =  1  in) 


Mold  Fraction 


•  Not  modeling  purge  could  partially 
explain  earlier  disagreements 
®  Axisymmetric  purge  computed  here 
is  only  partially  realistic 

■  Real  purge  was  on  only  one  side  of 
a  square  chamber 

■  Realistic  model  requires  3D 


HjMole  Fraction  (X  =  51n) 


O,  Mole  Fraction  (X  =  5  in) 
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Mole  Frsctie 


Mole  Fraction 


Effect  of  turbulence  model 

nm . ___ . - . 


©  The  effect  of  the  turbulence  model  is  not  small 
©  Should  use  a  model  appropriate  for  shear  layers 


,,  Summary  of  gas/gas  CFD 


•  Few  shortcuts  appear  to  be  possible  if  one  is  to 
truly  compare  experiments  with  CFD. 

©  Range  of  variation  suggests  that  steady  and 
transient  3D  computations  should  also  be 
investigated 
h  Underway  at  AFRL 


Other  available  data  sets 


•  Focus  of  the  2nd  International  Workshop  on  Rocket 
Combustion  Modeling,  25-27  March  2001, 
Lampoldshausen,  Germany. 

h  Data  set  RCM  -1 :  Round  liquid  nitrogen  jet  injected  into 
gaseous  nitrogen  (no  combustion)  at  3.97  and  5.98  MPa. 

h  Data  set  RCM  -  2:  Coaxial  LOX/H2  combustion  at  1  MPa. 
h  Data  set  RCM  -  3:  Coaxial  LOX/H2  combustion  at  6  MPa. 


•  Conference  results: 

a  Reasonable  ability  to  model  RCM  -  1 . 
m  Questionable  ability  to  model  RCM  -  2  and  RCM  -  3. 


Injection  at  supercritical  pressures 


The  problem 


•  Rocket  combustion  chambers  often  operate  at 
pressures  exceeding  the  critical  pressure  of  one 
or  more  propellant,  such  as  LOX. 

©  At  supercritical  pressures,  the  distinct  difference 
between  gas  and  liquid  phases  disappears. 
u  Conventional  “spray  combustion”  experience  no  longer 
applies. 

•  It  is  not  known  how  to  replace  conventional 
“spray  combustion”  models  in  engine  design 
codes. 

u  The  lack  of  understanding  leads  to  potentially  large 
engine  design  errors. 


Surface  tension  a  vanishes  at  supercritical  conditions. 
Conventional  atomization  and  breakup  parameters  become 
infinite ,  where  no  data  exists. 


Supercritical  atomization  and  breakup  regimes  are  largely 

unknown 


^  The  problem 

Other  factors  not  normally  considered  in 
conventional  spray  combustion 
©  Vanishing  surface  tension  and  enthalpy  of  vaporization. 

©  Equivalent  “gas"  and  “liquid”  phase  densities. 

©  Strongly  enhanced  solubility  of  one  species  (“gas”)  into 
another  (“liquid”). 

©  Reduced  “gas”  phase  diffusivity  (more  liquid-like). 

©  Large  property  excursions  near  the  critical  point 
a  Conductivity,  viscosity,  speed  of  sound,  specific  heats. 

©  Mixing  induced  critical  point  variations. 

»  Enhanced  gas  phase  unsteadiness. 

»  Potentially  different  kinetics  mechanisms. 


^  AFRL  facility 


»  Windowed  pressure  vessel 
operating  at  supercritical 
pressures. 

©  Cryogenic  fluid  capability 
(LOX,  LN2) 

9  Capability  to  produce  super¬ 
critical  droplets  and  jets. 

®  Shadowgraph,  Schlieren,  and 
Raman  visualization  of  con¬ 
centration  fields. 

•  Capability  to  drive  flows  with 
an  acoustic  driver 


Low  Reynolds  number  jet  results 


Low  Reynolds 
Number 
Cryogenic  Jets 


Shadowgraph  images  of 
liquid  nitrogen  jets  is¬ 
suing  into  a  pressurized 
chamber. 

Inj.  diam:  0.25  mm 
Re:  3350  -  4090 
LN2  Temp:  87K 
Chamb.  Temp:  292K 


7.4  mm 

t° 

9.5  mm  L 


15.0  mm 
to 

17.1mm 


22.6  mm  (  -  j  jj,  .Ct-  'h, 

10  t  L,  ‘'-vl:;,' 

24.7mm  f.tQLi 

t  \r m 


■# 


0.83  1.03 


Empirical  model 


O  N2  jet  into  N2  L/D=200  (’) 

♦  Cold  He  jet  into  N2;  L/D=200  (') 

M  Cold  N2  jet  into  He;  L/D=200  (*) 

x."  Brown  &  Roshko  (He/N2)lncompressrble  variable-density  Mixing  Layer 
A  02  jet  into  N2;  L/D  =  200  (“) 

- Incompressible  Variable-Density  (Papnmoschou&Roshko  theory) 

- Dimotakis  (theory) 

- Proposed  Model  (Chehrou di  et  al.) 

S '4  N2  jet  into  He;  L/D=200  (*) 

X  N2  into  (0.1  CO+0.9N2);  l/D:=20O  (*) 

A  N2  into  (0.5CO+0.5N2);  L/D=200  (*) 

X  N2  into  (0.5He+0.5N2)  (*) 

— “ — Proposed  Model  (Chehroudi  et  al) 


Gas-like  behavior  of  supercritical  jets 


®  It  has  quantitatively  shown  that  supercritical  jets 
behave  like  low  Mach  number  variable  density 
jets  in  several  respects 
h  They  have  the  same  appearance 
«*  They  spread  at  the  same  rate 
«  They  have  the  same  fractal  dimension 

•  Thus  there  is  a  reasonable  expectation  that 
“conventional”  turbulence  research  on  low  Mach 
number  variable  density  jets  will  apply  to 
supercritical  jets 


Jet  deformation  mechanism 


Bag  breakup  of  drops 
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Falling  drop  Jet  in  acoustic  field 


®  Drops  flatten  perpendicular 
to  the  flow  before  breaking 
up  because  the  higher 
velocities  around  the 

shoulders  cause  lower 

pressures  there  due  to  the 
Bernoulli  effect. 

*  Similar  pressure 

imbalances  are  set  up,  in 
the  mean ,  around  a  jet  (or 
drop)  in  an  acoustic  field. 

•  Therefore  the  jet  deforms 
in  the  mean. 


Main  Conclusions  from  Acoustics 
Observations 


» The  effect  of  acoustics: 
m  is  largest  near  the  critical  pressure, 
a  is  strong  at  subcritical  pressures, 
a  decreases  to  negligible  at  supercritical  pressures 
The  effect  of  acoustics  decreases  at  all  pressures 
as  jet  velocity  increases  (residence  time  in  the 
acoustic  field  decreases) 

Conjecture 

Being  in  a  supercritical  regime  may  reduce  the 
ability  to  couple  with  acoustic  instabilities. 


Transcritical  LN2  drops  in  room  temperature 
GN2 


1/16”  (1.6  mm) 


: 
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Representative  evolution  of  transcritical  drop  disintegration 


Messages  from  beyond? 


Super-Optical  natural  droplets? 


